Interaction-induced oscillations of the tunneling density of states 
in a non-quantizing magnetic field 
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We study tunneling into interacting disordered two-dimensional electron gas in a non-quantizing 
magnetic field, which does not cause the standard de Haas - van Alphen oscillations. Interaction 
induces a new type of oscillations in the tunneling density of states with the characteristic period 
of cyclotron quantum hu> c . 
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It is well-known that strong magnetic field modifies the 
single-particle density of states (DOS) of non-interacting 
electrons due to the Landau quantization. In a two- 
dimensional electron gas the quantization leads to a peak 
structure in the DOS, which is revealed in tunneling ex- 
periments as peaks in the dependence of the tunneling 
conductance on applied bias, see, e.g., Ref. Q. The form 
and the width of these peaks are determined || by the 
disorder. Experiments [|l],^[^] show also suppression of 
the conductance at zero bias. This suppression is a man- 
ifestation of interaction between electrons and has been 
studied theoretically both for disordered g and for clean 
H two-dimensional conductors in the quantizing mag- 
netic field. 

In a weak magnetic field the distance between the Lan- 
dau levels, %u>ci is smaller than their disorder-induced 
width. In such "classical" magnetic field the Landau 
quantization peaks in DOS disappear, while the afore- 
mentioned interaction correction to tunneling DOS sur- 
vives and in the limit of zero field evolves into a singu- 
lar at the Fermi energy negative logarithmical correction 
predicted originally by Altshuler, Aronov, and Lee Jjj. 

The question arises, whether or not in such "classi- 
cal", u) c t s <C f, limit the magnetic field influences the 
spatially averaged density of electron states measured 
in tunneling experiments (here r s is the electron quan- 
tum lifetime). For the non- interacting system, effect of 
the weak magnetic field on DOS is exponentially small 
0, oc exp[— 2n/(id c T s )}, and can be neglected. The goal 
of the present paper is to show that, to the contrary 
to the non-interacting case, in the interacting electron 
gas the "classical" magnetic field does produce a signif- 
icant effect on tunneling DOS. This effect becomes pro- 
nounced if the disorder potential is weak enough and 
smooth, with the correlation length much larger than the 
Fermi wave length and the amplitude much smaller than 
the Fermi energy. In such a potential, electrons expe- 
rience small-angle scattering, and their transport relax- 
ation time r tr , is much larger than r s . Thus there exists 
range of magnetic fields, in which Landau quantization 
is suppressed (lo c t s <C 1), while classical electron trajec- 
tories are strongly affected by the field (w c Ttr 3> 1). In 
this regime interaction correction to the tunneling DOS, 



Si>(e), is strongly enhanced with respect to the zero mag- 
netic field case. Furthermore, we will show that it ex- 
hibits peaks as a function of energy with the distance be- 
tween peaks equal to the cyclotron quantum, hu c . The 
further a peak is away from the Fermi level, the smaller 
and wider it is. The shape of the n-th peak in DOS, 
|e — nhuj c \ < hu) c /2, is given by: 



8v{e) h(uj c Tt r ) 2 1 fe — rihuj. 



where 
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energy e is measured from the Fermi level, and v — 
m/nh 2 is the free-electron density of states (m is the 
electron mass). The peaks overlap strongly for e > 
tku cy /Lu c T tI , and the oscillatory structure is washed out. 

Sensitivity of tunneling DOS to the classical magnetic 
field comes from the fact, that, as we will show, the 
interaction correction to tunneling DOS is associated 
with the self-crossing of classical electron trajectories. 
We denote the probability for an electron to complete 
a loop of self-crossing trajectory over time t as K(t). 
The interaction correction to DOS, 5v(e), turns out to 
be related to the Fourier transform of this probability, 
5v(e) oc K(e) = J °° dte~ let K(t). The strong enough, 
w c r tr 2> L magnetic field curves the electron trajecto- 
ries, significantly affects the return probability and, in 
turn, affects the tunneling DOS. 

For long time scales t ^> r tr , the function K(t) can 
be found from the diffusion equation. It gives K(t) oc 
(Dt)~ 1 for the two-dimensional case (D is the diffusion 
coefficient). The Fourier transform, K(e), is proportional 
to ln(e), which leads to a predicted by Altshuler, Aronov 
and Lee logarithmic correction to DOS at small ener- 
gies, e <C Ji/T ir , with the renormalized by the magnetic 
field diffusion coefficient Q. 

At short time scales, f < T tt , electrons move ballisti- 
cally along the cyclotron orbits. Provided that w c T tr ^> 1, 
during the time t electron may return to the initial point 
many times. Thus at these short time scales the magnetic 
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field significantly increases the return probability K(t). 
Multiple periodic returns of electron produce peaks in the 
probability Fourier transform K(e) at energies, which are 
multiples of the cyclotron quantum. Tunneling density of 
states oscillates with the same period, which is reflected 
by Eq. <§. 

Now we derive expression for the interaction correc- 
tion to DOS valid for arbitrary energies. We will put 
h = 1 in all intermediate formulas. In an ideal pure sam- 
ple the electron density does not depend on coordinates, 
n e (r) = no. Due to the sharp Fermi edge, scattering of 
the electron forming the Fermi sea on impurities results 
in an interference pattern in electron density. This pat- 
tern is commonly referred to as the Friedel oscillation 
||. In general, one can express density profile of nonin- 
teracting electrons in terms of the exact retarded Green 
function, Qf (r,r') = ^,*(r%(r)/(e - e, + iO) , of an 
electron in the random potential: 



n e (r) 



2 f° 

— / delm£f(r,r). 

J- ep . 



(3) 



Single electron wave function ?Az( r ) satisfies the 
Schrodinger equation for noninteracting electrons, 
H ipi = (e, + E F )in, where H = -(h 2 /2m)V 2 + U r (r), 
Ep is the Fermi energy, and U r (r) is the random poten- 
tial. In the presence of interaction V(r — r') between 
electrons, the Friedel oscillation produces an additional 
term in the Hamiltonian, Hhf, which can be presented 
(see e.g. Ref. as a sum of Hartree, Vh, and exchange, 
Vf, terms: 

H HF {r, r') = V H (r) S(r-r') + V F (r, r') (4) 
V h {t) = J V(r-r")Sn e (r")dr" (5) 



V F (r, r') = - V (r - r')Sp(r, r') 



(6) 



Here 5p{v, r') is the perturbation of the density matrix, 
p(r, r') = — (2/7r) J delmQf (r, r'), by the random poten- 
tial. Only the electrons with the same spin participate in 
the exchange interaction, which is reflected by the factor 
1/2 in Eq. (g). The Hartree-Fock energy (||)-(f|) oscillates 
as a function of coordinate in the same manner as Sn e (r) 
does. 

The local DOS is related to the retarded Green func- 
tion of electron, v{e,r) = — (2/tt) ImQ^(r, r). Scattering 
of electron on the Hartree-Fock potential, Eqs. 
induces a correction to Q^(r,r), which, in the Born ap- 
proximation, can be expressed as: 



^f(r,r) 



0f(r, r')T^(r')£f(r',r)dr' 



J gf{T,T')V F {v\v")gf{v",v)dv'dv". (7) 



We will be interested in the spatially averaged density of 
states, 8v{e) = (1/5) J 6v(e,r)dr, where S is the area of 



the system. For simplicity we will start with the case of 
the finite-range interaction potential and will calculate 
the Hartree contribution to the averaged DOS. Making 
use of Eqs. (0), (||), and (Q), and exploiting the identity, 
/ Qf(r', r) gf (r, r")dv = dgf{r', r")/de, we obtain: 



Su H (e) = J da J dr'dr"V{r' - r") 



dc 



(8) 



where gf(r,r') = [g^(r' ,r)]* . We are interested in the 
correction to the density of states averaged over the re- 
alizations of the disorder potential. Average of the prod- 
uct of two retarded Green functions g^(r', r')0*(r",r") 
does not contain contributions associated with the elec- 
tron trajectories longer than Xp. Thus, this product does 
not produce DOS energy dependence at small, as com- 
pared to Ep, energies, and can be neglected. On the 
contrary, averaged product Q R Q A is determined by long 
electron trajectories. Furthermore, it can be expressed 
JlOt in terms of the classical probability density T>: 

(^ 1 (r 1 ,r 2 )^ 2 (r 3 ,r 4 ))=^y —J — 

xel p 1 (r 1 -r 4 ) ei p 2 (r 3 -r 2 ) I3(ei _ ^ ^ ^.^ ^ ); 

V(u;l;2)= dte luJt V(t;l;2). (9) 
Jo 

Here p.; = pi?(cos0i, sin<^), and X>(t;l;2) is the proba- 
bility density for electron which starts at moment t = 
in point ri with the direction of momentum cf>i to arrive 
at moment t to the point r 2 with momentum direction 

4>2- 

Eq. (||) is valid as long as the sizes | r i — r 4 1 and | r 2 — r 3 1 
of spatial domains defining the ends of a trajectory are 
small enough, so that electron propagation in these two 
domains can be described by plane waves. Thus formula 
(||) is valid if e F r t ^ 3> 1 (semiclassical regime), and for 
the arguments sufficiently close to each other pairwise; 
in general, the condition |ri — r 4 |, |r2 — <C vpr s must 
be satisfied, however, in the special case of Eq. (ph, we 
need the product of Green functions with ri = r 3 and 
r 2 = r 4 , and this condition can be eased. Indeed, the in- 
tegral over the angular variables in Eq. (^) is dominated 
by the close to each other momenta pi and p 2 , and there- 
fore a weaker requirement, |ri — r 4 |, |r 2 — r 3 | <C hr,R c , 
should be satisfied. Here R c — vf/uj c is the cyclotron 
radius, l tI — v F r tI is the electron transport relaxation 
length, and v F and p F are the Fermi velocity and mo- 
mentum respectively. Further derivation of the Hartree 
correction requires substitution of Eq. (^), with argu- 
ments ri = r 2 = r' and r 3 = r 4 = r", into the Eq. (8) 
and integration over the difference between r' and r". 
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These coordinates are coupled by the interaction poten- 
tial V(r' — r"). Hence performing the integration we can 
exploit Eq. @ only if the range of the interaction poten- 
tial d is sufficiently short, d <C hr, R c - 

For a macroscopically homogeneous sample, classical 
probability T> depends only on difference of its coordi- 
nates r = Y2 — V\. Using this fact and Eqs. (||) and (^), 
we find: 

fag(e) _ 2 Re r°°^ f dfadfa dV(co; r = 0; fa, fa) 



x V(2p F |sin[(0i-0 2 )/2]|), 



(10) 



where V(q) is a Fourier transform of the interaction po- 
tential. Formula (|l^) can be easily generalized for the 
case of the long range Coulomb potential V(q) = 2ne 2 / nq 
(here k is the dielectric constant). For such a potential, 
V{q) in Eq. @ should be replaced (see e. g. Rcf. 0) 
by the screened potential I4 cr (q) = 2ne 2 /[n(q + 2/ag)] 
with as = h 2 K/me 2 being the effective Bohr radius. 
Note, that the range of screened Coulomb potential, 
V SCT (r), is of the order of the effective Bohr radius, and, 
therefore, much smaller than l tI and R c . Integration over 
the frequency in Eq. (|l0| ) immediately gives now the re- 
sulting expression for the Hartree contribution: 



5u H {e) 



-Re 

7T 



(2tt) 2 

x V SCT (2p F \sm[(fa 



V{e;r = 0;fa,fa) 
02)/2]|), 



(11) 



The exchange correction to DOS can be obtained in a 
similar way as Eq. ([ll]). The difference is, however, that 
in the exchange counterpart of Eq. (|l^), one should put 
the retarded screened potential V sct (uj, q) instead of V(q). 
The retardation of the interaction potential makes imme- 
diate integration over frequencies, that led to Eq. (11), 
no longer possible. The resulting expression for the ex- 
change correction has a form: 

8u F {e) 1 r [ dq , 0A(u;;q) 

- z Re/ dw/T^xaW^.q) — — , (12) 



V 7T 



<9w 



where F scr (u>,q) = V(q)/[1 + V(q)U(q,uj)}, with 
H(q, lo) — v[\ + iu;A(a;;q)] being the polarization op- 
erator, and 

A(o,;q) = /^/ ^dre-^ r ©(a;;r;^,0 2 ). (13) 

Equations ( |TT| ) and (|l^) express the interaction correc- 
tion to tunneling DOS in terms of classical probability 
density T>(e; r; fa, fa) and are valid at energies e % Ep. 
Function T>(e; r; fa, fa) can be found from the Boltzmann 
equation describing the scattering of electrons on impu- 
rities. In the special case r s <C r tr we are interested in, 
scattering on small angles dominates the collision inte- 
gral. With account for this simplification, the transport 
equation takes the Fokker-Planck form: 



— iio + 



P2 d 



d 1 d 2 



m e dr c dfa r tr dcj>\ 
2nS{fa-fa)6(r). 



V(uj;r; fa, fa) 



(14) 



Equation (|T^) describes electron motion along the cy- 
clotron orbit accompanied by angular diffusion caused by 
scattering on a random potential. Analysis of Eq. ( |l4] ) 
yields: 

f(w;q; 0i,<fo) = ^2'D n (uj;q;,fa,fa), (15) 



V n (uj;q; fa, fa) 



e m(^2-0l)gi-Hcq[(pi-P2)/PF xz] 

-i(w - nu> c ) + R 2 c q 2 /2T tr + n 2 jr tx ' 

(16) 



where z is a unit vector parallel to the magnetic field. The 



solution dl5[)-([l6D is valid for qR c < L0 2 T t i/(\ui\+uj c ). The 
obtained solution of transport equation together with 
Eqs. (|l|) and (|l^) enables us to calculate the inter- 
action correction to DOS in the classical magnetic field. 

At small frequencies, ui <C h/T tr , the n = term in 
Eq. (16) contains a part independent on initial, fa, and 
final, fa, directions of the electron momentum. This 
part dominates in Eq. (|i~5|), so that T>(oj; q; fa, fa) » 
l/(— iuj + R 2 q 2 / 2r tI ) ■ This limit corresponds to the diffu- 
sion regime studied in Refs. @,|]. Comparing Eqs. (|f) 
and ([ll]), one sees that in the diffusion regime the ex- 
change correction to DOS contains the interaction po- 
tential, V scr (uj, k), at very small, kR c < 1, momentum 
transfers, while for the Hartree correction the potential 
at large, k < p F , momentum transfers is important. 
Provided that the interaction potential decreases rapidly 
with the increase of k, one concludes JjJ that in the dif- 
fusion regime, e -C TtyV-tr, the exchange contribution to 
DOS, Eq. (|l~2|), dominates over the Hartree one. The re- 
sulting expression for the interaction correction to DOS 
in the this regime yields: 



8u(e) _ Ti(uj c T tT f 



87re F T tr 



In 



en. 



In 



\e\r tT a% 



TiR* 



(17) 



Equation ([l7]) corresponds to the Altshuler-Aronov-Lee 
M result with the renormalized by magnetic field diffu- 
sion coefficient D = R 2 /2rn- 

At larger frequencies, lu 3> h/T tl -, the probability 
T>(u>; q; fa, fa) describes quasiballistic motion of electron 
along the cyclotron orbit. After the period 2tt/oj c elec- 
tron approaches the vicinity of the initial point with a 
momentum only slightly deflected with respect to the ini- 
tial direction. According to Eq. ([ll]) , that means that in 
the quasiballistic limit the Hartree contribution to DOS 
contain interaction potential at small momentum trans- 
fers. Thus, in contrast to the diffusion regime, there is 
no special reason why the Hartree contribution should be 
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smaller than the exchange one. Furthermore, the Hartree 
correction appears to dominate over the exchange one. 
Indeed, when comparing the magnitudes of the Hartree 
and the exchange correction to DOS, one should take 
into account difference in a way how screening influences 
the two corrections. The Hartree correction contains 
screened interaction potential at zero frequency. On the 
other hand, expression for the exchange correction, Eq. 
(|l2|), contains integration over frequencies. The polariza- 
tion operator, H(uj,q), exhibits singularities at the same 
frequencies as the return probability, which causes sup- 
pression of the interaction potential V scl -(uj,q), and, in 
turn, suppression of the exchange contribution at ener- 
gies close to multiples of the cyclotron quantum. As a 
result, it is the Hartree contribution Jll| that determines 
the peak structure of the interaction correction to tun- 
neling DOS. Calculation of this correction, which con- 
sists in substituting Eqs. ([l5])-(|l6]) and expression for 
Vsci into Eq. ([ll]) and straightforward integration, gives 
the resulting expressions (|l|) and (||). Peak structure in 
tunneling DOS is well-pronounced for n < y?u) c T tY . For 
larger n width of the peaks becomes comparable with 
hu> c , and the oscillating structure of DOS disappears. 

Resulting energy dependence of the tunneling DOS of 
the interacting electron gas in a classical magnetic field, 
obtained by numerical integration of Eqs. ([n]) and fli~2] ) 
with account for Eq. (|l5|), is shown in Fig. 1. 
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FIG. 1. Energy dependence of the exchange (dashed), the 
Hartree (dotted) and the total (solid) interaction correction 
to tunneling DOS in a classical magnetic field. Energy e is 
measured from the Fermi level. The curves are calculated for 

l^cTtr = 40. 

The range of magnetic fields B where the oscillations 
of the density of states v(e) are caused by the interaction 
effects, is confined by the condition l/r tr <C oj c <C 1/t s , 
and depends on the sample quality. In a high mobility 
sample with the typical values ]l2] | of relaxation times, 
t s fa 10ps and r tr w 300ps this range is 1(T 3 T < B < 



5 • 10~ 2 T. Correspondingly, the typical numbers N of ob- 
servable peaks, N < yj rt r /r s , for these samples is N w 6. 

To conclude, we have shown that classical (uj c t s <c 1) 
magnetic field affects strongly the tunneling density of 
states of interacting electron gas. The tunneling DOS 
is found to be an oscillating function of energy with the 
characteristic period huj c . 
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